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ABSTRACT 



Context. Solar flares are powered by energy stored in the coronal magnetic field, a portion of which is released when the field 
reconfigures into a lower energy state. Investigation of sunspot magnetic field topology during flare activity is useful to improve our 
understanding of flaring processes. 

Aims. Here we investigate the deviation of the non-linear field configuration from that of the linear and potential configurations, and 
study the free energy available leading up to and after a flare. 

Methods. The evolution of the magnetic field in NOAA region 10953 was examined using data from Hinode/SOT-SP, over a period 
of 12 hours leading up to and after a GOES B 1 .0 flare. Previous work on this region found pre- and post-flare changes in photospheric 
vector magnetic field parameters of flux elements outside the primary sunspot. 3D geometry was thus investigated using potential, 
linear force-free, and non-linear force-free field extrapolations in order to fully understand the evolution of the field lines. 
Results. Traced field line geometrical and footpoint orientation differences show that the field does not completely relax to a fully 
potential or linear force-free state after the flare. Magnetic and free magnetic energies increase significantly ~ 6.5-2.5 hours before 
the flare by ~ 10" erg. After the flare, the non-linear force-free magnetic energy and free magnetic energies decrease but do not return 
to pre-flare 'quiet' values. 

Conclusions. The post-flare non-linear force-free field configuration is closer (but not equal) to that of the linear force-free field 
configuration than a potential one. However, the small degree of similarity suggests that partial Taylor relaxation has occurred over a 
time scale of ~ 3-4 hours. 

Key words. Sun: activity - flares - magnetic topology - photosphere - sunspots - corona 



1. Introduction 

Solar flares occur when energy stored in active region magnetic 
fields is suddenly released. The stored magnetic energy is con- 
verted to kinetic energy of particles, mass motions, and radi- 
ation emitted acro ss the entire electromagnetic spectrum (see 
iFletcher et al.ll20ri ). Energies of large flares can be as high as 
~ 10^^ ergs, over short time scales of tens of minutes. The field 
configuration of active regions is believed to be linked to the 
likelihood of flaring; flare triggers are often associated with flux 
emergen ce and increased s hear and twist in the field (see re- 
view by Priest & Forbes 2002, and references therein). However, 
the processes involved in magnetic energy storage and release 
within active regions are still not fully understood. Energy re- 
lease is expected to occur in the corona, but magnetic field ob- 
servations are usually based in the photosphere. In order to in- 
vestigate the coronal magnetic field, photospheric measurements 
are used as a starting point for 3D magnetic field extrapolations 
riarvlll989h . Although there are inaccuracies involved with the 
various assumptions that must be made to obtain 3D extrapo- 
lations, much progress has been made in recent years with high 
resolution photospheric magnetic field measurements now a vail- 
able from spacecraft such as Hinode (|Kosugi et alJl2007 h and 
Solar Dynamics Observatory jPesnell et alj 20121) . Increasingly 
accurate 3D magnetic field extrapolations have yielded insight 
into active region magn etic field topo l ogy during periods of flare 
activity (Regnier & C anfieldl 120061: iThalmann & WiegelmannI 
l2008tlSun et al..,2012) . 



This paper investigates the evolution of the 3D coronal mag- 
netic field in an active region using three types of extrapola- 
tion procedure: potential, linear force free (LFF), and non-linear 
force free (NLFF). The coron a is considered to be generally 
force free (iGold & Hovlell960l) . dominated by the relatively sta- 
ble magnetic field in a low-/3 plasma. Hence, all three types of 
3D extrapolation assume the force-free approximation, jx B = 0, 
where j is the current density and B is the magnetic field (iGarvl 
2001). This assumption is made under the special conditions of 
magnetohydrostatic equilibrium. The approximation results in 
the current being proportional to the magnetic field, with a pro- 
portionality factor a termed the force-free field parameter There 
are three generalised forms of the force-free relation, 



V X B = 0, 



V X B = aB, 



V X B = a(x,y,z)B. 



(1) 
(2) 
(3) 



A potential field configuration is defined as one containing no 
currents, resulting in the case of Eqn. [1] where a = 0. When 
a is non-zero but constant throughout a given volume the field 
configuration is referred to as LFF (Eqn. |2). Finally, when a is 
allowed to vary spatially (i.e., differing from field line to field 
line, but constant along one field line) the field configuration is 
referred to as NLFF (Eqn. |3). This specific case allows for the 
existence of both potential and non-potential fields within the 
given volume. 
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Potential field extrapolations can be thought of as the first- 
order approximation to the coronal magnetic topology. They are 
often used to represent solar global magnetic fields, which are 
dominated by sim ple dipolar configurations with few currents 
dLiu & LinI l2008h . In contrast, LFF field extrapolations have 
been considered sufficient to represent the large-scale, current- 
carrying coronal magnetic fields present on whole active re- 
gion size scales (Gary 1989; Wheatland 1999). However, NLFF 
field extrapolations intrinsically contain more information on the 
complex nature of solar magnetic fields and more accurately rep- 
resent coronal magnetic structures on size scales small e r than 
whole active regions. For example, Wiegelm ann et al.l (l2005l) 
compared field extrapolations from photospheric measurements 
to observed chromospheric magnetic fields, finding the potential 
solution provided no agreement with any observed field lines, 
the LFF solution provided 35% agreement, while the NLFF so- 
lution provided 64% agreement. It should be noted that there are 
still inconsistencies between existing NLFF extrapolation meth- 
ods, particularly regarding the trea tment of boundary co nditions 
(see, e.g., the comparison paper of iDe Rosa et al 

The underlying principle of force-free coronal fields ha s pre- 
viously been extensively studied. Initially, Wo ltiej (Il958l) theo- 
retically examined a magnetic field configuration over the course 
of magnetic energy release, proposing that force-free fields with 
constant a (LFF) represent the state of lowest magnetic en- 
ergy in a closed system. Tavlor (1986) applied this theory to 
laboratory plasma experiments, suggesting that the total mag- 
netic helicity of a flux system is invariant during the relaxation 
process to this minimum-energy state. The theory and experi- 
mental evidence led to the concept that the free magnetic en- 
ergy that may be released during field relaxation in a solar ac- 
tive region is the excess energy above the LFF fi eld with the 
same magnetic helicity (iHevvaerts & Pries3ll98 4V This relax- 
ation process will henceforth be referred to as Taylor relaxation 
throughout this paper Several theoretical, numerical, and obser- 
vational studies have investigated whether Taylor relaxation oc- 
curs during solar flares, with some agreei ng to its presence (e.g., 
iNandv et al.''2003'; 'Browning et al."2008") and others disagreeing 
(e.g.. .Amari & Luciani.2000: Blevbel et al.. 2002) . 

Investigating the magnetic and free magnetic energy in ac- 
tive regions is essential to study the physical processes occurring 
during solar flares. Previously, decreasing magnetic ener gy has 
been reported after flares (e.g.. lRegnier & Canfieldll2"006b . Also, 
energy budge ts for combined flare-CME events were studied in 
detail by Em sfie et al.1 (l2005h . finding -20-30% of the available 
free magnetic energy was required to power the events. 

In this paper, we examine active region magnetic energy evo- 
lution before and after a flare (on much shorter timescales than 
previously studied), in order to gain a better understanding of 
the processes involved in energy release. This paper compares 
the results of three different forms of 3D magnetic field extrapo- 
lation by examining the evolution of an active region. In Sect. |2] 
we briefly discuss the observations and analysis techniques em- 
ployed. Section m presents the main results, in particular geo- 
metrical differences in traced field-line solutions (Sect. 14.1b . ori- 
entation differences in field-line footpoint solutions (Sect. 14.2b . 
and the evolution of magnetic energy and free magnetic energy 
(Sect. 14.31 ). Finally, our conclusions and directions for future 
work are outlined in Sect. |5] 



2. Observations 

Active region NOAA 1095^ was observed by the Hinode Solar 
Optical Telescope (SOT: iTsuneta et al.l 1200 81 as it crossed the 
solar disk on 2007 April 29. The primary instrument used in this 
work from the SOT suite is the spectropolarimeter (SP). SOT- 
SP records the Stokes /, Q, U, and V polarization profiles of 
the Fei 6301.5 A and 6302.5 A lines simultaneously through a 
0.16"' X 164" slit. The data utilised here coiTespond to SOT-SP 
fast mapping mode (i.e., 0.32" x 0.32" pixels due to co-adding 
of adjacent slit positions and 2-pixel rebinning along the slit on- 
board the spacecraft). Details of the complete analysis proce- 
dure applied to the data are contained in Murray et al. (2012|) 
and references therein, including the atmospheric inversion pro- 
cedure, 180°-azimuth disambiguation, and heliographic coordi- 
nate conversion. Four SOT-SP scans are examined from 2007 
April 29 covering a period of 12 hours, leading up to and after 
a GOES Bl.O solar flare that peaked at 10:37 UT. Three scans 
were obtained before the flare (starting at 00:17 UT, 03:30 UT, 
and 08:00 UT, respectively) and one after the flare (starting at 
11:27 UT). Note that each scan takes 32 minutes to be built 
up from stepping the slit over the chosen field-of-view (FOV). 
Supplementary information on the location of the flare bright- 
ening is obtained from the Ca ii H filter of the SOT Broadband 
Filter Instrument. 

The active region investigated by [Murray et al.l is shown in 
the top row of Fig. [T] which is further investigated in this pa- 
per These panels contain the continuum intensity (upper left) 
and surface vertical field strength (upper rig ht) of the third scan 
(i.e., the scan prior to the flare). Previously. iMurrav et al.l iden- 
tified two regions of interest (ROIs; depicted in their Fig. 2) in 
an area of increased CaiiH flare emission. These regions are 
contained within the box in both upper panels of Fig. [T] (with 
green contours showing the CaiiH intensity at the flare peak). 
This zoomed-in box (also shown in the lower row of Fig. [TJ is 
used later for a portion of the analysis in Sect. 14.31 The surface 
magneti c field evolution of this region was previously consid- 
ered by Murray et al.l (l2012l) . finding significant changes in vec- 
tor field parameters in the ROIs during the observation period. 
Both ROIs showed a change in field inclination towards the ver- 
tical before the flare, with the field returning towards the hori- 
zontal afterwards. This variation in inclination agr ees with incli- 
nation changes across the neutral line predicted bv lHudson et al.l 
(2008). However, 3D coronal extrapolations are necessary to 
fully understand the evolution of the field, which is the specific 
purpose of this paper. 



3. Field Extrapolation Method 

In this work the output vector-magnetic field data analysed 
by iMurrav et al.l is subjected to further analysis. The helio- 
graphic vector field information (i.e., B'J, B^. and B^) obtained 
by Murray et al. were used as inputs to the three forms of 3D 
magnetic field extrapolation outlined in Sect. [T] It is necessary 
to preprocess photospheric magnetic field data before its use as 
a bounda ry condition for 3D extr apolations. The procedure out- 
lined by Wiegelm ann et al.l (12006) is used to deal with the in- 
consistency between the force-free assumption of NLFF models 
and the non force-free nature of the photospheric magnetic field, 
while also removing certain noise issues (such as uncertainties in 
the transverse field components). Note that preprocessing is ap- 



' http://solamionitor.org/region.php ?date=20070426&region= 10953 
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Fig. 1. Upper row shows the active region pre-flare state in continuum intensity (left) and vertical field strength (right). Green 
contours overlaid on both images indicate Ca ii H flare brightening observed at flare peak. The ~ 50" x 40" sub-region indicated by 
the black box in the upper row is shown in the lower row as preprocessed vertical field strength, at increasing scan time from left to 
right. A region of interest, thresholded at -750 G, is defined by the blue contours, and the purple line indicates the time of flaring 
between scans 3 and 4. 



plied to the data prior to the application of all three extrapolation 
methods. 

The pote ntial and LFF fields were calculated using 
themethod of Seehafer (1978), implemented in the LINFF code 
developed by T. Wiegelmann (2011, private communication). To 
obtain the potential extrapolation using LINFF, a is simply set 
to zero. For the case of LFF extrapolations, the required val- 
ues of co nstant a were calc ulated by fitting vs. from the 
results of iMurrav et al.l (see iHahn et al.ll2005h . The values of a 
obtained for the full FOV of each scan (i.e., 455 x 455 pixels^) 
are [0.30, 0.29, 0.27, 0.28] Mm"' for scans 1, 2, 3, and 4, respec- 
tively. The values are reasonably similar, showing a slight de- 
crease over time before the flare with a marginal increase after- 
wards. The chosen form of NLFF extrapolatio n is the weighted 
optimi zation method origin ally proposed by [Wheatland et aP 
(I2OOOI) and implemented bv IWiegehnannl (12004 ). This is cur- 
rently one of the most accurate NLFF procedures available, as 
den ionstrated in th e Schrijver et al. (2006), Metcalf et al. (2008), 
and lDe Rosa et alj (120091) method reviews. The NLFF optimiza- 
tion method begins by computing a potential field in the com- 



putational volume using the ISeehafej (Il978h method. The bot- 
tom boundary plane is defined by the input vector magnetogram, 
with the lateral and top boundary planes described using the po- 
tential field results. It is worth noting, this NLFF code directly 
minimises the force-balance equation, which avoids the explicit 
computation of a. 

A useful consistency check for the preprocessing method is 
the comparison of flux b alance, net force, and net torque param- 
eters (see Section 2 of Wiegelmann et al. 2006, for more details). 
In order to serve as a suitable bottom-boundary condition, vec- 
tor magnetograms must be approximately flux balanced, and net 
force and net torque must vanish relative to the force due to mag- 
netic pressure (Low 1985). The results for the data set used in 
this paper are promising. For example, the sum of net force and 
net torque normalised to magnetic pressure before preprocess- 
ing was on average ~ 0.4790, and after preprocessing this value 
drops to ~ 1.6 X 10"^. Thus, the preprocessing results in a data 
set more consistent with the assumption of a force-free coronal 
magnetic field. The bottom row of Fig. [T]shows the evolution of 
the preprocessed surface vertical field strength in the area of in- 
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Fig. 2. TRACE 195 A log-scale image of NOAA active region 
10953 at 03:45 UT on 29 April 2007. The overplotted white Hnes 
represent selected field lines from NLFF extrapolation results of 
scan 2. 



creased Ca ii H intensity (zoom-in region delineated by the box 
in the upper panels of Fig. [TJ. 

As mention ed in iDe Rosa et al ] d2009h and 

IWiegelmann et al ] (120 12h . it is useful to compare the ex- 
trapolated field with coronal-loop observations as a consistency 
check. This quantifies the extent to which the extrapolation 
correctly reproduces the coronal magnetic field configuration. 
Unfortunately very few high-resolution coronal observations 
were available for this particular event. Observations were 
obtained from the Transition Region and Coronal Explorer 
(TRACE: iHandv et alJil999) for comparison; a 195 A image 
was recorded at 03:45 UT on 29 April 2007 (i.e. , during scan 
2). Th e region of flare brightening examined by ' Murray et al.l 
(1201 2h contains only small-scale coronal loops, therefore 
large-scale loop structures from the full FOV were examined 
instead. Figure |2] shows select NLFF field lines traced over the 
image. Note that only the NLFF extrapolation is examined, as 
it is considered to be the most accurate representation of the 
coronal field (compared to the potential or LFF extrapolations). 
It is unfortunately difficult to identify coronal loops in this 
195 A wavelength image (171 A for example would have been 
more useful), but the extrapolated field lines seem to emanate 
from the leading sunspot to connect with areas of bright coronal 
emission. Agreement is particularly clear in the loops traced to 
the northeast of the image. 

A single ROI in the area of flare brightening was selected 
for further investigation. This was identified by thresholding the 
preprocessed surface vertical field strength at -750 G (depicted 
by contours in the lower panels of Fig. [T]i- Th i s RO I corre- 
sponds to 'ROI 2' investigated by iMurrav et al.l (1201 2h . but is 
not completely identical due to the data preprocessing (i.e., ef- 
fective smoo thing). This paper does not investigate 'ROI 1' of 
[Murray et al.l as it was deemed too small and fragmented to ob- 
serve any statistically meaningful variations in the extrapolated 



field. Extrapolation solutions were calculated for each of the 
four scans in computational volumes comprising of 455 x 455 x 
228 pixels^ (i.e., 105 x 105 x52 Mm^). It should be noted that the 
zoomed-in region (as shown in Fig. [TJ for further consideration 
in Sect. l4.3l makes use of the full height of the computational vol- 
ume, covering 165 x 135 x 228 pixels^^ (i.e., 38 x 31 x 52 Mm^). 
Also note that the values of LFF a obtained for the zoomed-in 
FOV of each scan are [0.23,0.26,0.22,0.17] Mm ' for scans 1, 
2, 3, and 4, respectively. 

4. Results 

Figure [3] shows zoom-ins on the results of the three types of ex- 
trapolation from the full FOV, with increasing scan time from 
left to right (the flare occurs between scans 3 and 4). Field-line 
traces for every 30* pixel in the ROI identified in the bottom 
panel of Fig. [T] are shown. The potential and LFF field-line so- 
lutions look similar, while the NLFF field lines reach greater 
heights and connect further South. 

In order to study relaxation of the coronal field, it is neces- 
sary to determine differences in the field structure between the 
three forms of extrapolation solution. As noted in Sect. [2] the 
NLFF code used here does not compute a. An estimate of a 
for all pixels within the computational volume can be obtained 
using Equation [3] (V x Z? - aB) with the results of the NLFF 
extrapolation procedure. 

Average values for a were calculated for each traced field 
line within the full computational volume as well as in the 
ROI, and histograms of the results are shown in Figure For 
NLFF extrapolations a is assumed constant along the entire 
field line, and this allows examining the traced field line av- 
erages from the source pixels within a region. A consistency 
check was carried out on full field lines traced from the source 
pixels, and a was found to be constant within + 0.02 Mm"^ 
Figure [4| shows a to be very similar for all scans in the full 
computational volume, with the distributions peaking between 
~ 0.03 - 0.09 Mm ', and a slight variation in scan 3 towards 
larger values. For the ROI, a also seems to be similar between 
scans but a larger variation is observed, and distributions peaking 
between ~ 0.03 - 0.13 Mm The distribution for the third pre- 
flare scan dififers the most from the other scans, being shifted to- 
wards larger values of a. This indicates an increase in the amount 
of twist in the field in the hours leading up to the flare. The fourth 
scan distribution indicates a decrease in a, and hence twist, after 
the flare has occurred. 

Although the changes in a observed in the full and ROI 
FOVs over the four scans are not large, it is promising to find 
expected increases (decreases) before (after) the flare. A num- 
ber of other 3D magnetic field vector parameters will be exam- 
ined here, showing more significant changes over the observa- 
tion period. In particular, geometrical differences will be inves- 
tigated between field-line traces that are caused by the differing 
presence of currents in each extrapolation solution. Section 14.11 
presents an analysis of 3D spatial offsets between the extrapola- 
tion solution field lines, while diff'erences in field-line footpoint 
locations are presented in Sect. 14.21 Finally, the total magnetic 
and possible free magnetic energies are calculated in Sect. 14.31 

4.1. Spatial Differences in Field-line Traces 

To obtain an indirect indication of how a may be varying in the 
NLFF extrapolation solution, a first step is to determine to what 
degree the extrapolation solutions differ when considering field- 
line traces from the same starting pixel. This can be quantified 
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Potential LFF NLFF 




Fig. 3. Zoomed-in FOV of extrapolated volume (see Fig. [TJ at increasing scan time from left to right. Upper to lower row: potential 
field extrapolation, LFF extrapolation, and NLFF extrapolation. Note that the FOV is the same as the ~ 50" x 40" box shown in 

Fig.m 



in terms of offsets in the computational x, y, and z spatial coor- 
dinates at specific points along the length of a field-line trace. It 
should be noted that field lines traced from the same source pixel 
(i.e., from the ROI) may have different path lengths in each ex- 
trapolation solution through 3D space before they return to the 
photospheric boundary. In order to calculate 3D displacements 
between comparable locations along each field line, the variable 
length arrays of (x, y, z) coordinates for each field line trace were 
interpolated to 11 points. These correspond to relative locations 
separated by 1/10 of the total path length. Displacements were 
calculated in each of the x, y, and z spatial coordinates between 
the same relative locations (i.e., 1/10 along the potential trace 
minus 1/10 along the LFF and NLFF traces, etc. . . ) with the val- 
ues along one entire field-line trace being averaged. The larger 
this displacement value is the further apart the LFF and NLFF 



traces are from the potential case, and the more current there is 
in the system. 

Figure |5] shows the results of differencing and averaging 
these arrays of x, y, and z interpolated coordinates for each ROI 
originating field-line trace. Histogram values for all four scans 
are overlaid for comparison. Note that all histograms throughout 
the paper have been normalised to the total number of pixels of 
the ROI, so as to show the percentage number of pixels. In terms 
of X coordinates, the mainly positive distributions in the upper 
left panel show the potential traces generally reach further East 
than the LFF traces (i.e, reaching smaller x pixel values). The 
distributions for the (NLFF - POT) and (NLFF - LFF) cases are 
more broad (due to the NLFF solution effectively having a dis- 
tribution of a values), with the NLFF field lines not reaching as 
far East as the LFF field lines. In terms of y coordinates, the LFF 
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Fig. 4, Distribution of calculated NLFF a values from aU pixels 
within the full computational volume (upper) and ROI (lower). 
A bin size of 0.02 Mm"' was used. Ordinate show percentage 
occurrences, while the abscissa axis indicates the value of a in 
Mm"' . Each scan distribution is coloured as per the legend. 



field generally reaches further South than the potential configu- 
ration (shown by the mainly negative y-coordinate distributions 
in the upper middle panel). The central and lower-middle panels 
again show broader distributions, with the NLFF field-line traces 
reaching much further South than the potential and LFF traces. 

The z-coordinate distributions for (LFF - POT) are much 
narrower than those of the (NLFF - LFF) and (NLFF - POT) 
distributions, with a greater percentage of values close to zero. 
The height differences between the NLFF and both the potential 
and LFF field-line traces indicate that the NLFF field lines exist 
at larger heights on average within the computational volume. 
The z-coordinate heights were investigated further, by calculat- 
ing the maximum field-line height for all pixels within the ROI 
(see Fig. |6]l. No significant changes are observed with time for 
the potential or LFF extrapolations. The largest changes are ob- 
served in the NLFF results, with the apex height increasing be- 
fore the flare in scan 3, and then decreasing afterwards in scan 4. 
Note that NLFF fields reaching greater heights has been previ- 
ously found (Regnier & Priest 2007b; Liu et al. 2011). The char- 
acteristics derived from Fig. |5]are summarised in Fig. Q which 
illustrates typical field-line traces from an ROI pixel where the 
NLFF field has a larger magnitude of a (and hence a greater de- 
gree of twist) than that of the LFF field. 

Considering the variation of Fig. |5] across time, none of the 
coordinate difference distributions vary significantly from scans 
1 to 3, whereas larger changes are observed in scan 4 (i.e., after 
the flare). This is clearest for the y coordinates when considering 
the central and lower-middle panels. A portion of each of these 
distributions have shifted towards smaller values, with a greater 
percentage of the (NLFF - LFF) distribution having values of 
zero compared to the (NLFF - POT) distribution. Thus, after the 
flare some of the NLFF field lines lie closer in 3D space to the 
LFF field lines than to the potential ones. This confirms that the 
LFF field contains current (as it deviates from the current-free 
potential field) but that the NLFF field has stronger currents (as 
it deviates further). Smaller differences between the NLFF and 
LFF fields in the post-flare scan indicates a decrease in NLFF 
currents relative to those in the LFF field. However, some of this 
is due to the increased LFF a in scan 4 (Sect. |2]i. Even though 
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Fig. 6. Distribution of apex heights from all ROI pixels for po- 
tential (upper), LFF (middle), and NLFF (lower) extrapolations. 
Y-axes show percentage occurrences, while the x-axis indicates 
the height in Mm. Each scan distribution is coloured as per the 
legend. 




Fig. 7. Cartoon of typical extrapolated field line traces. The black 
area indicates the negative polarity ROI, and white area the plage 
region. The grids outline single pixels. Field-line traces and po- 
sition angles are indicated in blue (potential), orange (LFF), and 
green (NLFF). 



portions of the extrapolation solutions are similar after the flare, 
the NLFF field contains stronger currents than the LFF field. 
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Fig, 5. Differences in x (left), y (centre) and z (right) coordinates of field-line traces between LFF and potential (upper), NLFF and 
potential (middle), and NLFF and LFF extrapolations (lower). Each scan distribution is coloured as per the legend, while bin sizes 
of 1 .0 were used for all x and y coordinates, 0.5 for (LFF - POT) z coordinates, and 0.75 for (NLFF - POT) and (NLFF - LFF) z 
coordinates. 



4.2. Orientation Differences Between Field-line Footpoints 

As indicated in the previous subsection, the amount of current 
present in the coronal magnetic field affects the direction of field- 
line traces and ultimately the location of their footpoints (result- 
ing from the twist that currents introduce in the magnetic field). 
In order to get a better picture of the effective distribution of a 
values in the NLFF solutions, the footpoint position angle, (p, of 
a field-line trace is defined as, 

^.to„-i(ZlZZi], (4) 

\X{ - Xi) 

where x and y are pixel coordinates on the solar surface, sub- 
script 'i' denotes initial coordinate, and subscript 'f denotes fi- 
nal coordinate. The initial coordinates are taken as each pixel 
within the negative polarity ROL such that the footpoint position 
angle measures the counter-clockwise angle from Solar West to 
the traced field-line end footpoint in the positive polarity plage 
(see Fig.Qfor an illustration). 

Figure |8] contains the results of this calculation for each of 
the extrapolation formats from all four scans. The position an- 
gle distribution for the potential case (0^°^; upper panel) does 
not vary much over the first three scans, but a greater portion of 
the distribution occurs at larger angles after the flare. The po- 
sition angle distribution for the LFF case (0^^^; middle panel) 
is very similar to 0^°^, but shifted to larger angles by ~ 10-1 5° 
in all four scans. This is expected with the introduction of twist 
in the field from the presence of currents that are equally dis- 
tributed throughout space for the LFF case. Finally, the position 
angle distributions for the NLFF case ((p^^^^; lower panel) are 



shifted further and broadened. This is also expected, with the 
NLFF case having a distribution of effective a values throughout 
space (and hence a distribution of currents and twist in the field). 
The three pre-flare scans exhibit very similar <p^^^ distributions. 
However, in scan 4 after the flare the <p^^^ distribution has de- 
creased occurrence of lower angle values and shows a stronger 
peak at a similar position angle to the peak of the ^''^^ distribu- 
tion. It should be noted that significant portions of the <p^^ and 
^NLFF distributions do not overlap. 

The footpoint orientations and temporal variation reported so 
far correspond to non-localised characteristics of the magnetic 
field. The spatial distribution of differences between footpoint 
orientations from the different extrapolation solutions contains 
information on what portions of the ROI harbour strong cur- 
rents (and hence magnetic energy). For each ROI source pixel, 
footpoint position angle values for the three extrapolation traced 
field-line solutions were subtracted from one another to give the 
footpoint position angle difference, A(p. This calculated quantity 
is presented in Fig.|9] where larger values correspond to greater 
deviation in footpoint position angle between the selected pair 
of extrapolation solutions (i.e., greater twist in the field). The 
A0^^^"^°^ values are generally small over all scans, with a slight 
increase in scan 4. This again confirms the known increase of the 
LFF constant a value in that scan (Sect.|2|i. Larger quantities are 
observed for A^A^lff-pot ^j^^j ^^nlff-lff^ ^j^j^ ^j^g i^^^^t lo- 
cated in the South East of the ROI (i.e., nearest to the magnetic 
neutral line with the positive plage). This indicates that the por- 
tion of the ROI containing the strongest currents (thus magnetic 
energy) is that closest to the neutral line involved in the flare. 
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Fig. 8. Distribution of footpoint position angles from all ROI 
pixels for potential (upper), LFF (middle), and NLFF (lower) 
extrapolations. Y-axes show percentage occurrences, while the 
X-axis gives the end footpoint position angle in degrees counter- 
clockwise from Solar West. Each scan distribution is coloured as 
per the legend, while bin sizes of 10° were used throughout. 



Figure [TO]represents the footpoint position angle differences 
as histograms to aid in a more quantitative comparison. The up- 
per row clearly shows the fairly constant, narrow distribution of 
Aip^^^^'^^ over scans 1 to 3, which broadens and shifts to larger 
difference values after the flare (as indicated in Fig. |9). Both the 
^^NLFF-POT (middle row) and A^^^^-i^^ (bottom row) distri- 
butions are considerably broader, showing some minor variation 
over the three pre-flare scans. However, both distributions have 
a greater percentage of ROI pixels with angle differences close 
to zero in scan 4 after the flare. The offset observed between 
the LFF and potential cases in Fig. |8]results in the /^^nlff-lff 
distribution being shifted to lower values by ~ 10-1 5° from the 
^^NLFF-POT distribution in all scans. Combining both of these 
results, the NLFF and LFF post-flare fields appear to be closer 
in configuration than the NLFF and LFF pre-flare fields were. 

4.3. Magnetic and Free Magnetic Energies 

It is important to investigate the amount of energy stored in a 
magnetic configuration, as some of this stored energy is released 
in driving the flaring process. The analysis up to now has been 
concerned with studying geometrical differences between field- 
line traces in the three extrapolation solutions. However, the ef- 
fects that the potential, LFF, and NLFF values of a have on the 
coronal energy content can be calculated directly from the ex- 
trapolat ion solutions. T he magnetic energy, E^, contained in the 
field is (ISchmiddll964 . 



where V is the 3D computational volume under consideration 
and B is the magnitude of the magnetic field vector at every point 
in computational (x, y, z) space. However, the values of for 
each of the extrapolation solutio ns does no t correspond to the 
total energy available for flaring (lAlvll 19841) . Instead, two forms 
of free magnetic energy, AE„i, can be calculated. 



» pNLFF-POT _ pNLFF _ pPOT 



^£.NLFF-LFF _ £.NLFF _ ^LFF 



(6) 
(7) 



r 



(5) 



which correspond to using either the potential (Eqn. |6) or the 
LFF (Eqn. |7]l field configuration as the minimum energy state. 
Note that there is no free magnetic energy in a potential field 
configuration, but the LFF and NLFF configurations have free 
energy due to the presence of currents. Values for AE]f^^^'^^ are 
not computed here since the ultimate aim of this work is to study 
the relaxation of the NLFF field over the course of a flare. 

Two different volumes are used for the energy calculations: 
the whole active region volume (455 x 455 x 228 pixels-'), and 
the surface zoom-in represented by the boxes in Figs. [T]and |3] 
(i.e., 155 X 135 x 228 pixels-'). Note that the zoomed-in volume 
uses the whole height of the extrapolation volume. The zoomed- 
in box was chosen as it corresponds to the location of greatest 
flare emission in Ca ii H. The box extent was determined by first 
choosing the furthest footpoint coordinates of all ROI field-line 
traces from the four scans as boundary points, but a number of 
small magnetic flux elements enter and leave this region over the 
course of the observations. In order to avoid any spurious effects 
these transient flux elements might have on the calculated energy 
evolution, the box was enlarged by -15 pixels on each side. This 
results in minimal amounts of flux entering or leaving the box 
over the four scans (Fig. [T] lower panels). Note, the zoomed-in 
region excludes the entire sunspot umbra and penumbra. 

The temporal variation of the calculated E^ values are dis- 
played in the upper panels of Fig. (TT] with results from the 
whole active region volume presented in the left column and 
the zoomed-in volume in the right column. Note that error bars 
in this Fig. were determined using estimated uncertainties of 
+ 10 G for the LOS field and + 10 G for the transver se field on 
the photospheric boundary (Wie geknann et al]l2012l) . The per- 
centage errors on the photospheric boundary were assumed to 
extend upwards, and percentage errors in Z?x, By, and Z?z, for all 
3D space were combined. Understandably, values of from 
the whole active region volume are much greater than those 
from the zoomed-in volume. In addition, Zij^j^^^ is consistently 
larger than E]^, which is consistently larger than E^^. This 
corroborates the findings of Sects. 14.1 l and |4~2| such that for all 
scans the NLFF extrapolations contain more twist in the field 
(i.e., stronger currents and greater magnetic energies) than the 
LFF extrapolations. For both volumes under consideration, E^ 
varies little between the first two scans, but increases in magni- 
tude by scan 3 (i.e., prior to the flare). After the flare, most of the 
magnetic energies decrease marginally. However, none return to 
the earlier pre-flare 'quiet' values. Differences in E^^^ between 
scans 3 and 4 contain changes in total magnetic energy due to 
the combined effect of the flare and evolution of the region over 
a ~3-4- hour period. Between these scans the whole active region 
volume shows a decrease in E^^^ of (9 + 2) x lO'" erg, with the 
zoomed-in volume decreasing by (5 + 1) x 10"'° erg. Both of these 
changes are likely to be dominated by the flare, given the essen- 
tially constant level of total potential magnetic energy observed 
over the flare. 

The results of the free magnetic energy calculations are 
presented in the lower panels of Fig. [TT] Similar to the total 
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Fig. 9. Spatial variation of footpoint position angle difference, Aip, from originating ROI pixel locations (see Fig. [10] for value 
distributions). Scan time increases from left to right, while rows depict values of A^''^^^''^^ (upper), A0^''^^"^°^ (middle), and 
^^NLFF-LFF (jower). A neutral line is plotted in grey for context, as seen in the lower row of Figure [1] 
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Distributions of footpoint position angle difference, A(f> (see Fig. |9]for spatial representation). Scan time increases from left 
while rows depict values of A0^^^"^°^ (upper; bin size 3°), A^*^^^^"™^ (middle; bin size 6°), and (lower; bin 



magnetic energies, no significant changes are observed for ei- 
ther volume between the first two scans. Both A£'JJj^^^"''°^ and 
AE'JJj^'^'^"^^^ increase from scan 2 to 3 (i.e., just prior to the 
flare). After the flare, AE'JJJ'"^^"^^^ decreases towards pre-flare 
'quiet' values, remaining slightly higher than the level observed 
in scans 1 and 2. A£'JJj''^^"^°^ also decreases, but not as much 
as AEjJj^^^"^^^. In terms of absolute changes over the flare, 
^^■NLFF-LFF (jecrcases by (1.6 ± 0.3) x 10^' erg in the whole 



active region volume, while the zoomed-in volume decreases by 
(7 + 2) X 10^*^ erg. Caution should be taken when interpreting 
changes in /^^nlff-lff ^ygj. ^^^^^ (e.g., it is known here from 
Sect. |2] that the value of LFF a is different before and after the 
flare). However, it is at least encouraging that the changes in free 
magnetic energy above the LFF state are within a factor of 2 of 
the absolute changes in total NLFF magnetic energy. This in- 
dicates that ~15-25% of the free magnetic energy that existed 
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Fig. 11. Magnetic energy (E^, upper) and free magnetic energy (AE^, lower) for each extrapolation type, as identified in the legends. 
Energies for the whole active region volume are shown on the left with those for the smaller zoomed-in volume shown on the right, 
while vertical bars denote eiTors. The purple vertical lines between the scans 3 and 4 mark the flare peak time. 



above the LFF state before the flare has been lost from the sys- 
tem over the course of the flare. 

5. Discussion and Conclusions 

The magnetic field configuration of a ROI within an active re- 
gion has been studied over the course of a GOES Bl.O mag- 
nitude flare. Geometrical differences between field-line traces 
through potential, LFF, and NLFF extrapolation solutions have 
been analysed, as well as their resulting magnetic energies. It is 
found that the general orientation of ROI field-line footpoints do 
not change significantly in the hours leading up to the flare, de- 
spite the ROI showing an increase in total magnetic energies. 
However, there are signatures of field redistribution after the 
flare that indicate incomplete Taylor relaxation: a portion (i.e., 
not all) of the NLFF field configuration becomes similar to that 
of the LFF field. Consideration of the magnetic and free mag- 
netic energies after the flare indicate that the region still has en- 
ergy available for further flaring. It is worth noting that a GOES 
B 1 .2 flare occurs in the same spatial region on 2007 April 29 at 
14:35 UT, less than 3 hours after the final SOT-SP scan analysed 
here. 

A number of previous works have debated the relevance of 
Taylor relaxation to the flaring process. IXmari & Lucianil (120001) 
use the presence of non-linearities in the post-relaxation state 
of numerical simulations to su ggest that Tay l or's t heory does 
not apply to flares and CMEs. IBlevbel et all (l2002i) reach the 
same conclusion using observations, finding that the relaxed 
post-eruption state is inconsistent with a LFF state. However, 
a large amount of helicity was eject e d from the region studied. 
In relation to this, iReenier & Priest! (l2007ah state that if helic- 
ity is not conserved (e.g., during a CME) then the minimum 
energy state can be a potential one. It is also not yet known 
whether helicity conservation is the only constraint on relaxation 
(Pontin et al. 2011). Although helicity is not calculated here, it 
should be noted that the studied flare event has no associated 
CME. 



In contrast, other studies have demonstrated the presence of 
Taylor relaxation. Numerical simu l ations of energy release in 
a coronal loop by iBrowning et al.l (l2008l) indicate that the re- 
laxed equilibrium s tate corresponds closely to a constant a field. 
iNandv et al.l 12003*) report the observational detection of a pro- 
cess akin to partial Taylor relaxation in flare-productive active 
regions that never achieve completel y LFF states within their 
observation periods. It is worth noting. iNandv et akl find that the 
relaxation process occurs on of the order of a week. The 12 hour 
period studied here is a considerably shorter time scale, which 
may contribute to only partial Taylor relaxation being observed. 
It has also been suggested that partial relaxation is perhaps to be 
expected from a magnetically complex system with flux emer- 
gence and cancellation to be accounted for (Pontin et al. 201 1). 

In the build up to the flare, a clear increase is observed 
in all magnetic energies and free magnetic energies, occurring 
~6.5-2.5 hours prior to the start time of the event. On aver- 
age, in the zoomed-in region volume the magnetic energy in- 
creases by ~ 7 X lO-"' erg and free magnetic energy increases by 
~ 2 X 10^' erg. In the whole active region volume, the magnetic 
energy increases by ~ 20x 10^** erg and free magnetic energy in- 
creases by ~ 15 X lO-'' erg. This may indicate a shorter time scale 
for flare energy input than has been previously observed (e.g., a 
gradual increase in magnetic energy over the course of a day 
before an M-class flare is reported by Thalmann & WiegelmannI 
'2008), but could be related to the low magnitude of the event 
studied here. It is worth n oting, in s tudying Hinode data asso- 
ciated with X-class flares, lJingetal.1 found no clear and 
consistent pre-flare pattern in the temporal variation of free mag- 
netic energy above the potential state. 

Considering changes over the flare, a marginal decrease is 
observed in most of the magnetic energies. Previous authors 
have been primarily concerned with reporting changes in the 
free magnetic energy above the potential state over solar flares. 
For example. Sun et al. (2012) find a decrease in ^^nlff-pot 
of ~3 X 10^' erg within 1 hour of an X2.2 flare (which they 
believe to be an underestimation), iThaknann & WiegelmannI 
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(l2008l) find a decrease of ~5 x 10^^ erg over an M6.1 flare, 
and lThalmann et al ] (120081) find a decrease of ~2 x 10^' erg af- 
ter a Cl.O flare (~40% of the available free magnetic energy). 
The general scaling of the flare event magnitudes agree with the 
Bl.O flare-related changes of (3 or 8) xlO^^ erg in A£Nlff-pot 
(for the zoomed-in and whole active region volumes co nsidered 
here, respectively). However, 'Regnier & PriestJ (l2007al) suggest 
that A£'j^''^^"^°^ gives an upper limit for the energy that can be 
released during large flares, while AEjJj^^^"^^^ is a good estimate 
of the energy available for small flares. In this work the value 
of A£'j;jL''''"^PP prior to the flare is ~(3 or 6) xlO^' erg (for the 
zoomed-in and whole active region volumes, respectively) with 
a corresponding decrease over the flare of ~(1 or 2) xlO^' erg. 
This indicates that ~20-30% of the free magnetic energy above 
the LFF state is removed during the course of the flare. 

In summary, it seems that the magnetic configuration of ac- 
tive region NO A A 10953 was did not fully relax to either a po- 
tential or LFF state after the Bl.O flare on 2007 April 29. In ad- 
dition, the energy budget remained sufficient to trigger another 
flare within 4 hours. Finding an active region in a partially re- 
laxed state after previous flaring may then be a good indicator 
of impending flare activity. Also, the increase of magnetic en- 
ergy on short time scales before a flare could be useful for near- 
realtime forecasting. However, the methods used to obtain these 
quantities are too computationally intensive to be currently ap- 
plied in near-realtime. 
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